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We use inelastic neutron scattering and molecular dynamics simulation to investigate the inter- 
play between the structure and the fast sodium ion diffusion in various sodium silicates. With 
increasing temperature and decreasing density the structure factors exhibit an emerging prepeak 
around 0.9 A" 1 . We show, that this prepeak has its origin in the formation of sodium rich channels 
in the static structure. The channels serve as preferential ion conducting pathways in the relative 
immobile Si-0 matrix. On cooling below the glass transition this intermediate range order is frozen 
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Alkali silicates are a paradigm for multicomponent 
glass forming systems in which even at typical melt tem- 
peratures the mobility of the alkali ions exceeds that of 
the Si-0 network by orders of magnitude 

una. m- 

though such a high mobility of the alkali ions should be 
reflected in the underlying structure, up to now, it is still 
unclear how the alkali atoms are built into the Si-0 net- 
work. A distribution where an increasing alkali content 
causes increasing homogenous disruptions of the Si04 
tetrahedal network of pure silica 
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does not agree with 
the observed drastic difference in mobility. In addition, 
this distribution is in conflict with the highly non-linear 
dependence of the viscosity on alkali concentration 0, 0] ■ 
These observations lead to ideas that propose the exis- 
tence ofpreferential ion conducting pathways in a Si-0 
matrix gQ. 

Here, we present evidence for the existence of sodium 
diffusion channels in the static structure. Using neutron 
scattering experiments and molecular dynamics simula- 
tions on various sodium silicate melts and glasses, we in- 
vestigate the structure that provides the mobility of the 
alkali ions. We demonstrate the accuracy of the molecu- 
lar dynamics model by our neutron data. Then we make 
use of the detailed information obtained by molecular 
dynamics to elucidate the observed intermediate range 
order as seen by neutron scattering. 

Sodium di- (NS2), tri- (NS3) and tetra- (NS4) sili- 
cate glasses were synthesized from ultrapure Na2C03 and 
Si02 powders by fusion at f500K. For the inelastic neu- 
tron scattering experiments glassy samples were encap- 
sulated in Pt cells giving an annular sample geometry of 
40 mm in height, 22.5 mm in diameter and a f .25 mm wall 
thickness. Measurements were performed at room tem- 
perature and on liquid samples at temperatures between 
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1200 K and 1600 K on the neutron time-of-flight spec- 
trometer IN 6 of the Institut Laue-Langevin in Grenoble. 
An incident neutron wavelength of A = 5.9 A -1 yielded 
an energy resolution of 5E = 50 fieV (FWHM) and an 
accessible momentum transfer at zero energy transfer of 
q = 0.2 — 1.75 A -1 . The raw data reduction consisted of 
a normalization to a vanadium standard and a correction 
for self absorption and container scattering. 

At 1600 K the structural relaxation of the Si-0 net- 
work for the investigated sodium silicates is on a nanosec- 
ond scale @,EJ - too slow to be resolved on IN 6. There- 
fore, the dynamic structure factor S(q,u>), depending on 
momentum q and energy transfer hu, exhibits a strong 
elastic line with an elastic structure factor S(q,uj = 0) 
that is in good approximation the static structure fac- 
tor times the Debye- Waller factor f(q). Whereas scat- 
tering on the Si and O atoms is exclusively coherent, 
scattering on sodium is coherent and incoherent. The in- 
coherent scattering reflects itself in a flat background in 
S(q,u = 0). 

The potential that was used for the molecular dynam- 
ics simulations is a slight modification of the pair po- 
tential by Kramer et al. |ll| which is based on ab ini- 
tio calculations. More details on this potential can be 
found in Ref. [l2|. The simulations were done at various 
densities for systems of N — 8064 particles (Nsi = 2016, 
7V Na = 1344, and N = 4704). At the densities 2.2g/cm 3 
and 2.37 g/cm 3 systems at the temperature T = 2100 K, 
i.e. in the liquid state, were equilibrated for 3.3 ns fol- 
lowed by microcanonical production runs of the same 
length (using the velocity form of the Verlet algorithm 
with a time step of 1.6 fs). In each case two indepen- 
dent runs were done in order to improve the statistics. 
The equilibrated systems at the density p = 2.37 g/cm 3 
were also quenched to a glass state at T = 300 K using a 
cooling rate of about 10 12 K/s. In the latter simulation 
the pressure was kept constant at ambient pressure which 
yielded a density of 2.46 g/cm 3 at 300 K. Note that the 
densities p = 2.2 g/cm 3 and p = 2.46 g/cm 3 are close to 
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the experimental densities at T = 2100 K and T = 300 K, 
respectively @. 

The central quantity for the determination of the elas- 
tic structure factor from the simulation is the time- 
dependent dynamic structure factor defined by 

S ( q > *) = V AT h2 Yl b ^ b 0Sa0(q, t) a, fie {Si, Na, 0} 

(1) 

with 

-, N a Nf, 

Sa/3(g,t) = ^£X;<e X p[iq-(r fc (t)-r,(0))]> , (2) 

k=l 1 = 1 

Yi(t) being the position vector of particle i at time t. In 
Eq. JiQ b a denotes the experimental coherent scattering 
lengths of particle species a. The elastic structure factors 
were estimated by the equation S(q,uj = 0) = f(q)S(q) 
(with S(q) = S(q,t = 0)) where the Debye-Waller fac- 
tor f(q) was determined by fitting in the fi— relaxation 
regime the normalized quantity F(q,t) = S(q,t)/S(q) to 
a von Schweidler law (for details see Ref. [ljj)- Note that 
we have not taken into account the incoherent scattering 
on sodium in the calculation of S{q,u> = 0) because this 
contribution is not expected to change the shape of the 
function S(q, u) = 0) significantly. 

Fig. displays the elastic structure factor of glassy 
and viscous sodium trisilicate as seen by inelastic neu- 
tron scattering. The maximum at ~1.7A -1 corresponds 
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to the disrupted tetrahedal Si-0 network structure of al- 
kali silicate 0, 0. With increasing temperature the 
elastic scattered intensity at this wave- vector is decreas- 
ing at the expense of an increasing inelastic scattering. 
In contrast S(q,Lu = 0) displays an additional prepeak 
at ~ 0.9 A" 1 which grows with increasing temperature, 
indicating an enhanced probability for a correlation of 
atomic arrangements at distances around 6-8 A. 

In sodium disilicate and sodium tetrasilicate neu- 
tron scattering reveals a similar behaviour: The elastic 
structure factors exhibit a growing prepeak at ~ 0.9 A -1 
on temperature increase. With increasing sodium con- 
tent the height of the prepeak is increasing, whereas its 
position remains fairly unaffected by this large change in 
composition (Fig. QJ)). We note, that for a monoatomic 
system a change in the structure that is due to thermal 
expansion would lead - beside a slight shift of the peak 
positions towards a smaller q value - to a decreasing in- 
tensity of the peak height with increasing temperature. 
Thus, the emerging prepeak in the elastic structure fac- 
tors as shown in Fig.^does not seem to be merely caused 
by thermal expansion but might be related to a change 
in the underlying structure. In the following, we clarify 
this behaviour of the elastic structure factors by means 
of our molecular dynamics simulations. 

Fig.|3displays the "simulated" neutron scattering elas- 
tic structure factor of sodium trisilicate. At the ex- 
perimental densities, S(q,u> — 0) exhibits the same 
behavior as in the neutron scattering experiment |16| : 
Whereas around 1.7 A -1 the elastic scattered intensity 
decreases with increasing temperature (as expected), 
around 0.9 A^ 1 a shoulder is present at T = 2100 K which 
is nearly absent in the glass at T = 300 K. Also shown 
in Fig. |2 is the elastic structure factor at the density of 
2.37 g/cm 3 at T = 2100 K. Here the shoulder around 
0.9 A -1 has a smaller amplitude than in the correspond- 
ing case at the experimental density. This gives evidence 
that the possible structural changes leading to a more 
pronounced appearance of the feature around 0.9 A -1 in 
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FIG. 1: Elastic structure factor as seen by inelastic neutron 
scattering: Note the emerging prepeak at ~ 0.9 A -1 with in- 
creasing temperature in sodium trisilicate (a) with a position 
that is fairly independent of sodium concentration (b). 



FIG. 2: Elastic structure factor of sodium trisilicate as ob- 
tained by molecular dynamics simulation (see text) weighted 
with the neutron scattering lengths: Note the emerging pre- 
peak at ~ 0.9 A -1 with decreasing density. 
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FIG. 3: Partial structure factors So,o{q), SWa,Na(5), Ssi,Na(g) 
in glassy and liquid sodium trisilicate at experimental densi- 
ties. So,o(q) + 0.1 for clarity. 

S(q, u> = 0) are not due to a change in temperature but 
due to the change in the density. Thus, the question 
arises whether the shoulder around 0.9 A -1 is related to 
structural features that tend to disappear with increasing 
density. 

The partial static structure factors S a p(q) at the ex- 
perimental densities for different correlations are shown 
in Fig. |21 A well-pronounced peak is present in <Ssi,Na(<z) 
and SWa.Na^) at T = 2100K and at T = 300K. Thus, the 
structure that leads to the peak at q\ = 0.9 A -1 does not 
at all disappear with increasing density. That the feature 
at qi seems to be absent in S(q, ui = 0) at 300 K, is due to 
the fact that S(q, u> — 0) is a linear combination of six dif- 
ferent partial structure factors: On the one hand positive 
and negative contributions cancel each other (note the 
negative amplitude of S'si,Na(o) at gi), and on the other 
hand So,o(q) exhibits only a shoulder at qi which is less 
pronounced at T = 300 K. Because oxygen is the major- 
ity component (about 60% of the particles in NS3) and 
the coherent scattering length of oxygen is significantly 
larger than that of sodium and silicon (6o/&Si ~ 1-4 and 
fro/^Na ~ 1.6), So,o(q) gives the major contribution to 
S(q,w = 0) and thus, the changes in So,o(q) are the 
main cause for the emerging prepeak in the experimental 
and simulated S(q,uj = 0) with increasing temperature 
and decreasing density, respectively. 

The feature at qi corresponds to a distribution of 
sodium ions that is not homogenous on a length scale 
of 6-8 A. This is illustrated by the snapshot, Fig. 0] 
which shows a molecular dynamics configuration of NS3 
at T = 2100 K at the density 2.2 g/cm 3 . Here, the sodium 
ions are represented by blue spheres that are connected 
to each other by surfaces. For clarity, the Si and O atoms 
are shown as small yellow and red spheres, respectively, 
that are connected to each other by covalent bonds shown 
as sticks. Note that the size of the latter spheres does not 
mirror the real size of Si and O. The snapshot gives a clear 



picture of a network of sodium channels that percolate 
through the Si-0 structure. The peak at qi in the static 
structure factor marks the characteristic length scale of 
the latter network of channels. 

Recently, Jund et al. p"?| have shown in a molecular 
dynamics simulation of NS4 that the sodium trajectories 
form a network of pockets and channels where the dis- 
tance between the pockets is of the order of 5-8 A. Then, 
Horbach et al. [18| have given evidence in the case of 
NS2 that the latter network is reflected by the prepeak 
at qi in the static structure factor. Fig.0]shows that the 
"structure of the sodium trajectories" is not only seen in 
the dynamics as claimed by Jund et al. • The sodium 
trajectories are given by paths in a quasi-static struc- 
ture of sodium diffusion channels, the lifetime of which 
is given by the characteristic relaxation time of the Si-0 
network. 

A change of the Na 2 content (from NS4 to NS2) or of 
the density does not significantly affect the position of the 
prepeak at qi as seen by molecular dynamics simulations 
and neutron scattering. The effect that the prepeak at qi 
becomes more pronounced in S^q, lo — 0) with increasing 
Na 2 concentration (Fig. QJj) can be also explained by 
the simulation [T^ : The amplitude of the peak around 
1.7 A -1 decreases with an increasing sodium content in- 
dicating a stronger disruption of the Si-0 network if more 
sodium ions are added to the system. Moreover, Na-Na 
correlations contribute with a larger weight in the total 
structure factor and thus also the amplitude of the pre- 
peak at qi increases. These findings from the simulation 
explain the behavior of S(q,uj — 0) as seen in Fig. ^i. 

However, a change in the composition or the density is 
accompanied by a significant change in the disrupted Si- 
O network structure. This is indicated by the behaviour 
of the Si-O-Si angle: In the NS3 simulation at the ex- 
perimental densities this angle increases from 134.5° at 




FIG. 4: Molecular dynamics snapshot of the structure of 
sodium trisilicate at 2100 K at the density p = 2.2 g/cm 3 : 
The blue spheres that are connected to each other represent 
the Na atoms. The Si-O network is drawn by yellow (Si) and 
red (O) spheres that are connected to each other by covalent 
bonds shown as sticks between Si and O spheres. 
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T = 300 K (i.e. at the density of 2.46 g/cm 3 ) to 143° 
(i.e. at the density of 2.2 g/cm 3 ). Thus, the sodium dif- 
fusion channels are sustained over a broad range in tem- 
perature and sodium concentration because the packing 
of the disrupted Si-0 tetrahedra is rather "flexible" with 
respect to the Si-O-Si angle. Of course, there should be 
a lower bound of sodium concentrations below which the 
latter rearrangements in the Si-0 network are no longer 
possible. Indeed, at very low sodium oxide content (be- 
low about 5mol%), a miscibility gap of Na 2 and Si0 2 
even at typical melt temperatures has been reported . 

Very recently, Lammert et al. |20| have investigated 
a lithium silicate melt by molecular dynamics simula- 
tion. They have shown that the diffusion mechanism for 
lithium is a kind of vacancy diffusion. Movies that we 
have done from the NS3 configurations such as Fig. 0] 
reveal a similar sodium diffusion mechanism. Thus, our 
picture of diffusion channels is consistent with the va- 
cancy diffusion on a local length scale as proposed by 
Lammert et al. provided that also in lithium silicates 
diffusion channels exist. Indeed, our recent results of in- 
elastic neutron scattering measurements |2l) and molec- 
ular dynamics simulations |22| on various lithium and 
potassium bearing silicates indicate that preferential ion 
conduction pathways are a more general feature for al- 
kali silicates. A neutron scattering experiment on a cal- 
cium silicate glass using Ca isotope substitution revealed 
a prepeak in the Ca-Ca partial structure factor [23j at 
~ 1.3 A -1 . In contrast to our findings in alkali silicates, 



the signal has been interpreted in terms of a glass struc- 
ture based on ordered, densily packed domains. It re- 
mains to be seen, to what extent the existence of prefer- 
ential ion conducting pathways in the static structure 
shown here for sodium silicates applies to other glass 
forming, oxide ion conductors. 

In conclusion, inelastic neutron scattering experiments 
and molecular dynamics simulations on sodium silicate 
melts and glasses revealed the existence of sodium dif- 
fusion channels in the static structure. We find, that 
the formation of preferential ion conducting pathways 
is fairly unaffected by the sodium oxide concentration 
(in the studied composition range from sodium tetrasili- 
cate to sodium disilicate), by a temperature change cor- 
responding to the glass at 300 K up to the equilibrium 
melt, or a change of the density in the melt of about 8 %. 
In sodium silicates the values of the sodium ion diffu- 
sion at typical melt temperatures 0, 0, Q an( A the values 
of the viscosities @, Q depend relatively weakly on the 
sodium oxide concentration. The formation of sodium 
rich channels gives an explanation on a microscopic level 
for the observed macroscopic properties of mass trans- 
port in sodium silicates. 
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